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The main concepts of the recently developed approach to singular problems of 
quantum mechanics are extended to the Dirac particle in the Coulomb field of a 
point-like nucleus with its charge Z >137. The reflection and transmission coef- 
ficients, which describe, respectively, the reflection of electron by the singularity 
and its falling onto it, are analytically calculated, using the exact solutions of 
the Dirac equation. The superradiation phenomenon is found. It suppresses the 
widely discussed effect of spontaneous electron- positron pair creation. 



1 Introduction 



In our previous publications [TJ E] we developed a new approach - called 
black-hole approach - to the Schrodinger equation with attractive inverse- 
radius- squared singularity in its potential, possessing the fall -down onto 
the center behavior. Within this approach the solutions, oscillating near the 
point of singularity, are treated as free particles, emitted or absorbed by the 
singularity exactly in the same sense, as they are emitted or absorbed by the 
infinitely remote point - we mean, income or outcome from/to infinity. 

In the present work we apply this approach to the relativistic problem of 
an electron placed in a Coulomb field of a point-like nucleus with its charge Z 
greater than 137. The corresponding set of radial Dirac equations is known to 
possess in this case the singularity of the same mathematical character as in 
the Schrodinger problem above. Namely, the Coulomb potential is squared in 
the Dirac equation and produces singular attraction (Za/r) 2 . Basing on the 
same concept and using the known explicit solution of the Dirac equation, we 
are able to calculate the coefficients of reflection and transmission for electrons 
that income from infinity and are then either reflected back or transmitted to 
the singularity point. In the negative energy continuum e < —m, where m is 
electron mass, these coefficients are responsible for the spontaneous electron- 
positron pair creation by the external field of the nucleus. This effect has been 
studied in a more conventional context (the basic references are [3], [3]), using 
an introduction of the nucleus finite-size core, acting as a cutoff. In the present 
approach the reflection and transmission coefficients are found analytically in 
a simple form, since they are not affected by a cutoff, whose presence is not 
meant at all. The result shows the superradiation phenomenon, characteristic 
of the Kerr black holes, when the reflected wave is more intense than the 
incident one. The reflection coefficient is greater, than unity, the transmission 
coefficient being, correspondingly, negative. This fact forbids the spontaneous 
pair creation owing to the Pauli principle applied within the Dirac sea picture. 

The paper is organized as follows. In Section 2 we shortly outline the 
concepts of the black-hole approach, addressing the Dirac equation in the 
supercritical Coulomb field, with the reference to the original work |3] for 
details. In Subsection 2.1 we formulate the Kamke generalized eigenvalue 
problem with respect to Z in what is called sector III of confined states 
— m < e < m, Z > 137; find the corresponding integration measure in the 
orthogonality relations, that allows treating the solutions, oscillating near the 



singularity, as free particles; perform the transformation (when < e < m) 
that reduces the Kamke eigenvalue problem to the standard Liouville form 
and find the corresponding effective potential, confining the particle near the 
singularity. In Subsection 2.2 we analogously formulate the effective barrier 
problem in what is called sector IV, e 2 > m 2 , Z > 137, where particles emit- 
ted by the singularity or by the infinitely remote point are reflected back or 
transmitted to, resp., the infinitely remote point or the singularity. In Section 
3 the reflection/transmission coefficients are calculated. The recipe for this 
calculation is determined by the concept, but is independent of any of the 
constructions in the previous Section 2. 



2 Outline of the black-hole approach 

The radial Dirac equation in the Coulomb field is two-component [U], [7j 

dG(r) k , . Za. . 

— + -G(r) -(e + m + )F r = 0, 

dr r r 

T v ' + (e-mH )G(r) F(r) = 0. (1) 

dr r r 

Here the spinor components of the radial wave function G(r) and F(r) cor- 
respond to gr and /r in notations of Ref.jE], resp., e and m are the electron 
energy and mass, a = 1/137 is the fine structure constant, Z is the charge of 
the nucleus, and k is the electron orbital momentum 

k = -(Z + 1) for j = l + ~, 

K = l for 3 = l~\- (2) 

In what follows we confine ourselves to the lowest orbital state k — — 1. Solu- 
tions to equation ([TJ are known in terms of confluent hypergeometric functions. 
To illustrate the connection of the problem (JIJ with the Schrodinger problem, 
which possesses the fall-down onto the center property, let us for example 
consider the combination 

It satisfies the second-order Schrodinger-like differential equation 



d 2 W fVm 2 -e 2 1 2 2 \ (Z 2 a 2 -\ 2Zae\ Tjr , A . 

+ ~ m )W=( -= +——)W. (4) 



We collected all terms containing the coupling Za in the right-hand side. 
This fixes the operator, whose eigenf unctions will be studied as solutions to 
the generalized Kamke eigenvalue problem jH] with respect to Za. The term 
Za/r in (j3J) corresponds to Coulomb attraction of an electron by a positively 
charged nucleus, when Za and e are positive, while the term Z 2 a 2 /r 2 produces 
singularity, which is attractive irrespective of the sign of Za. This property 
of the relativistic problem is well-known. The fall-down onto the center jU] 
occurs, provided that 

Z 2 a 2 - 1 > 0. (5) 
Two fundamental solutions to equation (0J) behave near the origin r — ► 

like 



where 7 = \/ Z 2 a 2 — 1. 

Two fundamental solutions to equation (J1J) behave near the infinity r —>■ 00 

like 



c 



TiZae ± 1 



The usual Coulomb bound states of an electron lie within the strip 

— m < e < m, aZ < 1, (8) 

corresponding to what was called sector I in 012112] • In this sector we choose 
the solution, which decreases exponentially at infinity (the upper sign in (JJJ)). 
It is a linear combination of two solutions, out of which one decreases near 
r = (7 is imaginary in sector I), whereas the other increases in accordance 
with (JHJ) and should be ruled out by the requirement of square-integrability. 
The nullification of the corresponding coefficient is possible for discrete values 
of e, considered as functions of aZ: e nr (aZ), n r = 0, 1, 2... 
On the contrary, in the strip 

— m < e < m, aZ > 1, (9) 

corresponding to what was called sector III in 0121 E|, the both behaviors © 
oscillate, and, in accord with the black- hole approach of Refs. jT] , [2] , |Sj we 
treat them as particles, free near the origin. These particles are emitted or 
absorbed by the singular center, but do not escape to infinity due to Q, with 



the upper sign again. The corresponding states make a continuum, because 
the both oscillating solutions are acceptable. We call them confined states in 
order to distinguish from the bound states that also do not extend to infinity, 
but belong to discrete spectrum. The treatment of the solutions, oscillating 
near r = 0, as responsible for free particles, is possible, since the measure in 
the corresponding scalar product becomes infinite in this point. Within the 
present context this will be seen later. 
In sector II 



electron and positron are elastically scattered from the Coulomb potential. 
Finally, in sector IV 



the wave function oscillates both near the origin and at infinity, the particles 
being free in the both regions. This corresponds to any of the two inelastic 
scattering processes: An electron (or a positron) incoming from infinity, is 
partially scattered back with the probability, determined by the reflection 
coefficient, and partially penetrates to the origin, becoming a free particle 
there. Symmetrically, the particle, emitted by the center, escapes to infinity 
with a certain probability, determined by the transmission coefficient, and 
partially is reflected back to the center. We say that deconfinement processes 
go in sector IV. 

2.1 Sector III of confined states — m < e < m, Za > 1 

Following [U 121 E] , equation (JTJ) or (jlj) should be in sector III considered as a 
generalized eigenvalue problem of Kamke [8] for the parameter Z, while its 
solutions should be classified as eigen-functions of the operator in the left-hand 
side of (HJ). Following the standard procedure, we may multiply eq.(jH) by W*, 
the solution to the equation, complex conjugate to (jU, written for a different 
eigen- value Z\. After subtracting the complex conjugate of the same product 
with Zi and Z interchanged, and integrating by parts, we obtain 



e < — m or e > m, 



aZ < 1 



(10) 



e < — m or e > m, 



aZ > 1 



(11) 




(12) 



where 




(13) 



is the integration measure with which two eigenf unctions with different eigen- 
values Z\ 7^ Z and the same e are orthogonal, provided that the boundary 
term on the left-hand side in (jl2j) disappears. This is guaranteed if the system 
is placed in a box 



tl < r < oo, ri 
and the appropriate boundary condition 

W(r L ) = 



0. 



(14) 



(15) 



is imposed at the wall of the box. 

Let us make the transformation of the variable ( ro is a free dimensional 
parameter) 



£(r) 



2Zae\ 



1 



dr 



(16) 



and of the wave function 

W(£) = $(r)W(r), where $(r) 



7 



+ 



2Zas\ 3 
r / 



(17) 



This transformation is nonsingular when the energy e > is positive, and 
reduces equation (J1J) to the standard Liouville form 



d 2 W(Q 



+ U(t)W(£) = - 1)W(0 



with the effective potential 



2(1 + 



2Zae r \2 



e 1 + r 2 (e 2 - m 2 ) + | 



(19) 



It is understood that r is here the function of £, inverse to (|16)). The free 
parameter r enters the potential only through this function. Unlike the non- 
relativistic case of pure inverse square potential (V = in Refs. |U 121 El), 
where there is no dimensional parameter originally present in the theory, so 
that it is to be additionally introduced to define the problem, in the present 
case under consideration we have at our disposal the Compton length m _1 . 

For negative e the potential (jTUjl becomes singular as a result of the singu- 
larity of the transformation (|16|). (|T7j) . 



In the vicinity of the singular point r = the transformation (|17j) 
behaves as 

C = ln^, (20) 

(I)V(r) = #(0 (21) 

The origin r = is mapped onto £ = — oo, and the infinitely remote point 
r = oo to the point £ = oo (£ x {2\r2Zae/^)^/r). The asymptotic behavior 
© in the point r = becomes the oscillating asymptotic behavior at £ — > — oo 
(up to an unessential constant factor) 

^(0xexp(±i£ 7 ) (22) 

with 7 = \/ a 2 Z 2 — 1 playing the role of a "momentum". The potential |T9|) 
grows at £ — > oo as ((m 2 — e 2 )j 2 /2aZe)£ 2 ( see Fig.l ) and prevents the 
electron from escaping to infinity. 
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Figure 1: The confining potential in sector III plotted for Za = 2. The 
lower curve corresponds to the positive energy value e = 0.5m, and the upper 
one to £ = 0.1m. The value r = m -1 is taken for r 



The box (fTf]) is mapped to the box , restricted by the limits 

- £ L < £ < oo, (23) 



where 



r L = r exp (-£l), £l -> oo, when r L -> 0. (24) 

The boundary condition (|T5)l becomes 

W(-Cl) = (25) 

in full analogy with the customary way of quantization of particles, free at 
infinity. Eq. (|23j) makes the Kamke problem self-adjoint and provides the 
orthogonality and completeness of the wave functions (see j3] for detail.) 

When necessary, the finite quantity may be considered as a cutoff, al- 
ternative to the nucleus core introduced in [I]. 

As long as the box size is finite, the spectrum of Za is discrete in the 
domain Za > 1 and is presented by an equidistant series of vertical lines, 
numbered by n = 0, 1,2... from left to right in the plane with y = e for the 
ordinate, and x = Za > for the abscissa, which fill the domain, restricted by 
the borders of sector III, i.e. —m < y < m, 1 < x < oo. As the spectral lines 
cannot intersect, these vertical lines must be adjusted to the Coulomb bound 
states levels, which are described in the same plane by a family of trajectories, 
numbered by the radial quantum number n r < N. N is finite, provided that 
tl is finite. These trajectories go in sector I, as Za grows, from the point 
y = m, x = to the border x = 1 , separating the two sectors. The lowest 
and leftmost trajectory n r = joins the leftmost vertical line n = near the 
point y = 0, x = 1 and is continued by it down toy = —m. Other trajectories 
with n r = 1,2...N join with the vertical lines n = 1,2. ..AT above that point 
also to be continued down to e — —m. The other spectral lines with n > N 
remain unaffected and go almost vertically from y = m to y = —m. When 
tl — ► 0, all the spectral lines to the right of the border x = 1 condense to make 
the continuum of confined states. In this limit the bound state trajectories 
approach that continuum at the border of sector III at x — Za = 1 and at 
discrete positions on the vertical axis, corresponding to the discrete energy 
eigenvalues, numbered by n r . 

When £ L — > oo, the norm of solutions of eq. (fTH|) 

/ \W(0\ 2 ^ = 1 ! \W(r)\ 2 % (26) 

diverges linearly with £x- This makes the ^-normalization of the wave functions 
possible. Note the singularity of the volume element dr/r 2 in the initial r- 
space: the particle has enough volume near the origin to behave itself freely 
near it. 

Such is the situation with sector III. 



2.2 Sector IV of inelastic scattering a 2 Z 2 > 1, e 2 > m 2 . 
Effective barrier 

Excluding F(r) from equation ((TJ, one comes j3] to the equation for the upper 
component of the Dirac spinor G(r) 



d 2 G Za dG (IZae -«(e + m) \ 

dr 2 r(Za + r(m + e)) dr \ r r(Z« + r(m + e)) / 

+ £ 2 - m 2 ) G. (27) 



We placed the most singular term in the r.-h. side, as well the "kinetic energy" 
term e 2 — m 2 . 

The corresponding equation for the lower component F(r) is obtained from 
this by the substitution G —>■ F, Z — > —Z, n — ► — «, e — > — e. 

Set k = — 1 and perform the transformation of the wave function G, pre- 
scribed by the general theory of differential equations jH], 

G({) = *(r)G(r), 
t(r)=f a ' Z l- 1 +^-"''Vf r( ^ +n 'l m + £ !lV. P8) 



r 2 / \ro(o;Z + r(TO + e)) 

accompanied by the transformation of the variable 

£(r) = / _ v / a 2 Z 2_ 1+r 2( £ 2_ m 2^ (29) 

</r r 

When r — > the new variable tends to — oo as lnr. When r — > oo the new 
variable tends to oo as £ — > (aA 2 — m 2 /^f)r. After the transformation (J2B1), 
equation (|27|) acquires the standard Liouville form 



+ t/(OG(0 = (ZV - 1)G(0 (30) 



d£ 2 

with the effective potential 



U(0 



— - 2aZer + — ^ + — — ^ ^ 



(l + y 2 )V aZ + r(m + e) 4 (aZ + r(m + e)) 2 

2/ = r 7- ( 31 ) 



4(1 +V 2 ) 2 ) 1 * a 2 Z 2 -l 



If the transformation of the coordinate is not used - only the wave function is 
transformed - equation (J3U|) . (JBTj) is reduced to eqs. (2,11), (2,12), (2,13) of 
Ref.g]. 

Equation (j3Uj) should be supplemented by, e.g., periodic or antiperiodic 
boundary conditions in a box — £ L < £ < £jy, £ L — > oo, £{/ — > oo, where 
£i is defined as (J2U) (c/. [3J). The nonrelativistic case, considered in [3], 
teaches us that, in sector IV, the complete and orthogonal sets of solutions 
of the generalized Kamke eigenvalue problem make families, for each of which 
a certain relation between the two parameters Za and e is fixed, while the 
eigenfunctions within each family are labelled by the remaining free parameter. 
This relation for the present relativistic case remains unknown to us. 




-10 -5 5 10 15 20 25 

% 

Figure 2: The barrier potential (|3*T|) in the upper continuum. The curve is 
drawn for the positive energy e = 2m and the charge twice the critical value 
aZ = 2. The value ro = m _1 is taken for ro 



In the upper continuum (a part of sector IV) e > m the potential U(£) (|3*T|) 
is plotted against £ in Fig. 2. The most interesting situation appears in the 
lower continuum e < —m. The potential U(£) (j3~Tj) is plotted for this case in 
Fig. 3. It has the barrier character and a singularity in the point £ S i ng = £(r sing ), 



where 



sing 



—aZ 

m + e 



> 0. 



(32) 



This singularity first appeared in equation (|27j) , resulted from the exclusion of 
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Figure 3: The barrier potential (J3TJ) in the lower continuum. The curve is 
drawn for the negative energy e = —2m and the charge twice the critical value 
aZ = 2. The value r = m~ x is taken for r 



lower spinor component F from equation due to the zero in the coefficient 
in front of F in (Q). Vice versa, one might exclude the upper spinor component 
G from (Q). Then one whould encounter the second-order differential equation 
like (|27j) with the coefficients singular in another point r' sing = (Za/(m — e), 
provided that e < m. 

The singularity in the point also appears in the part of the potential 
U 2) "responsible for the spin interaction", eq.(2,13) of Ref.flJ, although is not 
discussed there. We admit, that the superradiation phenomenon discussed in 
the next Section may be attributed to this singularity. 



3 Superradiation and the pair creation by the vacuum 



The contents of this Subsection are not directly dependent on the procedures 
developed above, since here we deal with known solutions, not with equations. 

In sector IV all solutions are meaningful, neither should be ruled out. Let 
us consider the two-component solution to equations (0), whose behavior near 
the origin contains one oscillating exponent 



G(r) x (2pr) [V ^ Z 
F{r) x (2pr) iV ^ Z 



'i + Va 2 ^ 2 - 

j imaZ 

V 

(^r + Va 2 Z 2 - 1 

j \mctZ 

V 



— 1 



£ \ 2 

1 + - , 

mJ 



m 



(33) 



where p = y/e 2 — m 2 . The same solution behaves at infinity r — > oo as 



G(r) x £e- ipr (2pr)- i ^ (— + lV + ^(2^)^ f— + 1 



F(r) x -i£e" ipr (2pr)" 



1 



?7i 



1 1 

5 + i^e^(2pr) ia f s (— - lV (34) 



The other independent solution differs from this one by the change of sign in 
the exponent in (Correspondingly the constant coefficients £ and Q in 

become different.) The solution, described by the asymptotic equations 
(jnnj), (pSt is a wave, falling from the infinity and reflected back to infinity, with 
a transmitted wave, absorbed by the center. The coefficient Q is responsible 
for transmission/absorbtion, and £ for reflection. As a whole, in sector IV, the 
inelastic process is described by two scattering phases and one mixing angle. 
These are arranged in a unitary 2x2 ^-matrix, the reflection and transmission 
coefficients being expressed in terms of these angles [TJ |21 El • 

To find the absorbtion and reflection coefficients, note that the conservation 
divj = of the current j = ip^fip in the Dirac equation implies that the flux 



i / (jn)r 2 dfi = FG* - F*G, 



(35) 



where n = r/r and Q is the stereoscopic angle, be independent of r. Eq. (j35j) 
is the Wronsky determinant for the set (JI}. By equalizing the two expressions 
for the flux (|35|) obtained near r = and r = oo, we obtain the current 



conservation property in the form of the identity 



\£\ 



+ 



2Va 2 Z 2 - 1 



ISI 2 ( 70=? - vW^T 



(36) 



Here the second term in the r.-h. side is the transmission coefficient, and the 
first term the reflection coefficient. Using the known explicit solution, whose 
asymptotic behavior is (}3*3*j) . we find for the latter coefficient: 



where 
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Figure 4: The reflection coefficient in the upper continuum (for positive en- 
ergies e > m). The curve is drawn for the charge twice the critical value 
aZ = 2. 



The reflection coefficient (|37j) is plotted against the energy in Figs. 4, 5. In 
the positive kinetic energy domain e > m the barrier, as shown in Fig.2, is neg- 
ative (attractive). Correspondingly, the reflection coefficient (|37|) is very small 



Figure 5: The reflection coefficient in the lower continuum (for negative en- 
ergies e < —m). The curve is drawn for the charge twice the critical value 
aZ = 2. 



(see Fig. 4). Its maximum achieved at the border e = m makes exp(— 27r7). 
The Coulomb center acts as a strong absorber: it reflects next to nothing 
and absorbs almost everything. On the contrary, in the negative energy range 
e < —m there is a repulsing barrier singular in the point (|52*j) (|3ip. as shown in 
Fig. 3. Correspondingly, the reflection coefficient is greater than unity, as it fol- 
lows from (J37|) and is seen in Fig. 5. The phenomenon of R being greater than 
unity is called superradiation. It is known for scattering off the Kerr (rotating) 
black holes ^U] and is associated with the fact that there appears singularity 
within the definition domain of the differential equation. It is referred to as 
manifesting the possibility to pump out energy from a black hole. 

In our present context, the superradiation forbids the electron-positron 
creation from the vacuum in the field of a point-like nucleus with Za > 1. 

To explain this statement, let us first remind jU |5] that this process may 
be traditionally understood within the second-nonquantized theory as follow- 
ing. The nucleus is supplied with a core, which acts as a cutoff. When, for 
sufficiently large charge Z, the electron level sinks into the negative continuum 
e < —m, an electron taken from the filled Dirac sea can perform the under- 
barrier transition to become a tightly bound state near the nucleus. At the 



same time, the hole in the sea left by that electron, becomes a positron. In 
this way, the pair, comprised of a free positron and a deeply bound electron is 
produced from the vacuum by the strong cut-off Coulomb field. The probabil- 
ity of this process is determined by the probability of the penetration through 
the barrier. 

Our consideration of the nonregularized, point-like nucleus modifies this 
idea in only one respect. For us, the electron, after it passes through the 
barrier, becomes a confined state, free near the singularity. The process is de- 
scribed by the transition coefficient T = 1 — R. Now, once the superradiation 
takes place, the number of electrons reflected back into the Dirac sea, exceeds 
the number of those taken from the sea. This is forbidden by the Pauli princi- 
ple, since all vacancies in the sea are occupied. In this way the superradiation 
supports the vacuum stability. 

4 Conclusion 

The mathematical basis for application of the present black-hole approach to 
the relativistic problem of an electron in the overcritical Coulomb field is not 
as firmly established as it was in the case of the Schrodinger equation 
Lacking is a discription of families of sets of complete orthogonal solutions to 
the Kamke eigenvalue problem in sector IV, where inelastic processes occur. 
Nonetheless, the adjustment of the black-hole approach into relativistic con- 
text undertaken in Section 2 demonstrates its conceptual relevance, sufficient 
for calculating, in Section 3, the reflection and transmission coefficients, re- 
sponsible for spontaneous pair creation. The establishing of superradiation in 
electron scattering off the Coulomb center with Z > 137, which suppresses 
the pair creation is independent of the way the Kamke problem is posed and 
solved. 
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